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Abstract 



Let Ti(d) be the space of complex hermitian matrices of size d x d and let H+(d) C TL(d) be the 
cone of positive semidefinite matrices. A linear operator $ : Ti(d\) — > TL{d2) is said to be positive 
if $[i? + (di)] C H + (d2). The concurrence C($; •) of a positive operator $ : TL(di) — > Jiidi) is a 
real- valued function on the cone H+(di), defined as the largest convex function which coincides 



with 2 \J O2 2 {<&{£,(,*)) on all rank 1 matrices ££* € H + {d 1 ). Here erf : H(d) R denotes the second 

symmetric function, defined by erf (^4) = X^<j MiMji where fii, . . . , fid are the eigenvalues of A. The 
concurrence of a bipartite density matrix X is defined as the concurrence C(&;X) with $ being 
the partial trace. 

A analogous concept can be considered for Lorentz-positive maps. Let L n C R n be the n- 
dimensional Lorentz cone. Then a linear map T : K m —> R n is called Lorentz-positive if T[L m ] C 
L n . For this class of maps we are able to compute the concurrence explicitly. 

This allows us to obtain formulae for the concurrence of positive operators having Ti(2) as 
input space and consequently of bipartite density matrices of rank 2. Namely, let $ : 7i(2) — * 
H(d2) be a positive operator, and let Ai, ...,A4 be the generalized eigenvalues of the pencil 
<72 2 ($(X)) — AdetX, in decreasing order. Then the concurrence is given by the expression 



C($;X) = 2\J a 2 2 (<E>(X)) — A2 det X. As an application, we compute the concurrences of the 
density matrices of all graphs with 2 edges. 

Similar results apply for a function which we call /-fidelity, with the second largest generalized 
eigenvalue A2 replaced by the smallest generalized eigenvalue A4. 

1 Introduction 

The concurrence is a scalar function initially introduced to quantify the entanglement of bipartite den- 
sity matrices describing the mixed states of 2-qubits [BJ . In [BJ , an explicit formula for the concurrence 
of 2 <g> 2 bipartite density matrices of rank 2 was obtained. In a subsequent paper [T3j , Wootters gener- 
alized this formula to 2 ® 2 bipartite density matrices of arbitrary rank. Further generalizations were 
achieved by Uhlmann pj]. He considered real- valued functions /(££*) on the set of pure states and 
introduced the convex roof of /, which is the largest convex extension f(p) to the set of all density 
matrices. Similarly, the concave roof is the smallest concave extension of /. Uhlmann derived an 
explicit formula for the convex roof of the function /(££*) = |£*0£;|, where O is an arbitrary anti-linear 
hermitian operator acting on the state vector £. He called this convex roof 0-concurrence. It then 
turns out that the concurrence for 2 ® 2 bipartite density matrices just equals the 0-concurrence for 
a special anti- linear hermitian operator acting on C 4 , and Wootters formula is a special case of 
Uhlmanns formula for 0-concurrences. Uhlmann derived a similar formula for the ©-fidelity, which he 
defined as the concave roof of the function /(££*) = |£*0£l- 

Rungta et al. [9j defined the /-concurrence of arbitrary bipartite density matrices as convex roof 
of the function /(££*) = 2\/<J2 (^1 (££*)), where 02 is the second symmetric function of a matrix 
and tr\ is the partial trace with respect to the first of the two subsystems. Since the /-concurrence 
is the unique natural generalization of the concurrence as defined in [6J, we will henceforth call it 
simply concurrence. Rungta and Caves [10] computed the concurrence explicitly for d ® d bipartite 
density matrices of isotropic states, i.e. convex combinations of the maximally mixed state and a 
maximally entangled state. Osborne [7] obtained a formula for the tangle, i.e. the convex roof of the 
function /(££*) = 402(^1 (££*)), of rank two density matrices, using essentially the main idea in [6j. 
Uhlmann [12] then went up one abstraction level and replaced the partial trace tr\ in the formula 
= 2y/&2(tri(£ti*)) by an arbitrary positive operator i.e. he defined the concurrence C(&;p) 
of a state p with respect to the operator $ as the convex roof of the function /(££*) = 2 y/<72 ($(£,£,*))■ 
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He then showed that when $ is a completely positive map of rank and length two, the concurrence 
C(<£>; •) can be reduced to the 0-concurrence for a suitable anti-linear hermitian operator O. 

As noted above, the 0-fidelity differs from the 0-concurrence by the substitution of the convex roof 
by the concave roof. By analogy, we are tempted to introduce a function called /-fidelity by replacing 
the convex by a concave roof in the definition of the /-concurrence. For a bipartite density matrix, the 
/-fidelity will then be the smallest concave extension of the function /(££*) = 2 •^/^(tri (££*)) defined 
on the pure states. The /-fidelity F($>; •) of a positive operator $ is then the concave roof of the 
function /(ft*) = 2^2(*(ft*))- 

In the present contribution, we generalize the ideas of [B] in a completely different direction. Our 

point of departure is that the second symmetric function is a quadratic form with signature (H 

— •••—). Such forms are intimately linked with the second-order cones, or Lorentz cones, which are 
defined in real vector spaces R" of any dimension. Like the cones of positive semidefinite matrices, 
the Lorentz cones are so-called self-scaled cones, i.e. cones of squares for some Jordan algebra pQ. The 
structure defined on R™ by the Jordan algebra permits to introduce notions like eigenvalues, rank, 
trace and determinant for arbitrary elements of R™. Once these notions are adopted, the Lorentz cone 
naturally appears as the cone of positive semidefinite elements and the second symmetric function 
as determinant. For this reason, the Lorentz cones are particularly well adapted to the study of 
concurrence. In this paper we never use the Jordan algebra explicitly, so the reader is not required to 
be familiar with this concept. 

Initially the concurrence is defined only for density matrices, but by homogeneity it can be extended 
to the whole cone of positive semidefinite matrices [12]. Once the determinant is defined on R", we 
can generalize the notion of concurrence in the following manner. Let K C E be a convex cone defined 
in some real vector space E, and let T be the set of its generators, i.e. points which lie on extreme 
rays of K. Consider a linear map iff : E — ► R™ such that iff[K] C L n . We call such maps K-to- 
L n positive. Then we can define the concurrence C(iff;-) : K — > R of W as the convex roof of the 
function 7(7) = 2-y/dct(W(7)). Here / is defined on T and det is the determinant in K™ defined with 
respect to the Jordan structure associated with the Lorentz cone L n . We will show that this definition 
generalizes the notion of concurrence as defined in [12] in the following sense. For any positive map 
$ there exists n £ N+ and a map iff sharing the input space E with $ and having output space R™ 
such that C(<I>; •) = C(iff; •). The role of the cone K in the definition of the concurrence of \& is played 
by the cone of positive semidefinite matrices. Similar definitions and results can be obtained for the 
/-fidelity. 

If the input space E of iff is R m and the cone K is the Lorentz cone L m , then we are able to compute 
the concurrence C(iff; •) and the /-fidelity F(if>; •) explicitly. Since the Lorentz cone L4 is isomorphic 
to the cone of positive semidefinite complex hermitian 2x2 matrices, this allows us to compute the 
concurrence and the /-fidelity of a positive map $ whenever its input space is the space H(2) of complex 
hermitian 2x2 matrices. This, in turn, yields explicit formulae for the concurrence and the /-fidelity 
of bipartite density matrices of rank two. As an application, we compute the concurrences and the 
/-fidelities of the density matrices of all graphs with 2 edges, as defined in [2]. 

The paper is structured as follows. In the next section, we provide the necessary definitions. We 
recall the definitions of concurrence as in [9] and [12] and provide similar definitions for the /-fidelity. 
We introduce the Lorentz cones and define the necessary functions related to their Jordan structure. 
Then we generalize the notions of concurrence and /-fidelity to relate them to this structure. In section 
3 we investigate the relation between bipartite matrices and completely positive maps and show that 
the concurrence and /-fidelity of positive maps can be reduced to the concurrence and /-fidelity of some 
K-Xo-L n positive map. Sections 4 and 5 contain the main theorems of the paper, namely the formulae 
for the concurrence and the /-fidelity in the case when the input space of the K-to-L n positive map 
is R m and K is the Lorentz cone L m . In section 6 we concretize these results to the case of positive 
maps with input space 7i(2) and bipartite matrices of rank two. In the next section we apply these 
results to the density matrices of all graphs with 2 edges. Finally, we summarize our results and draw 
some conclusions in the last section. 

2 Definitions and preliminaries 

For some vector space E, let ids be the identity operator on E. Denote by /„ the n x n identity 
matrix, by i the imaginary unit of the complex numbers, by int W the interior of a set W and by dW 
its boundary. 

For an n x n matrix A, denote by 02 (-<4) its second symmetric function Xa<fc<Kn ^k^i, where 
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Ai, . . . , A„ are the eigenvalues of A. The second symmetric function can be written as 



(To 



l<k=£l<n \ \k=l / k=l J 



Let Ti{d) be the space of complex hermitian matrices of size d x d and let H + (d) C Ti.(d) be the 
cone of positive semidcfinite matrices. The space TL{d) has d 2 real dimensions. We equip the second 
symmetric function on Ti.(d) with an additional upper index d to indicate the size of its input matrices. 
If A E H.(d), then tr(A 2 ) = (A, A) is the squared Frobenius norm of A. Hence the second symmetric 
function on Ti{d) becomes 

4(A) = \^(I d ,A) 2 -|>f fe ,,4) 2 j , (1) 

where Mi, . . . , M c p is any orthonormal basis of the space Tt(d). 

A linear operator $ : H(di) — > W(d 2 ) is said to be positive if &{H + (di)] C H + (d2). It is said to 
be completely positive if for any n E N+, the operator id«( n ) £g>$ : H(n) ® H(dx) — > H(n) (g> ^(cfe) is 
positive. The maximal rank achieved by matrices in the image of <& is called the rank of <f>. 

A completely positive operator can always be represented as a sum 

N 

$(x) = j2 A « XA *k> 

k=l 

where the Kraus operators A% , . . . , An are complex di x d\ matrices. The minimum number N necessary 
for such a representation of <f> is called the length of $ [12] . 

Fix two numbers d\ , d-2. E N+ . A d\ ® d% bipartite matrix M is a matrix in 7i{d\d2) equipped with a 
block structure, namely partitioned into d\ x d\ blocks Mm of size di x di each, where the indices k, I run 
through 1, . . . , d%. Positive semidefmite d\ (g) d2 bipartite matrices with trace 1 describe the mixed state 
of a composite quantum system consisting of a subsystem having d± states and a subsystem having c?2 
states. The partial trace tr\ of a d\ <S> c?2 bipartite matrix M with respect to the first subsystem is the 
e?2 x di matrix Y^k=i Mkk- The partial trace tr2 with respect to the second subsystem is the d% x di 
matrix having tr Mm as (fc,/)-entry. 

The following definition of concurrence for positive operators is from [12j . 

Definition 2.1. The concurrence C(<i>; •) of a positive operator $ : TC(di) — » TL(d2) is a real-valued 



function on the cone H + (di), defined as the largest convex function which coincides with 2y er^ 2 ($(X)) 
on all rank 1 matrices X E H+{d\). 

We now introduce the following similar notion. 

Definition 2.2. The I-Edclity F(Q; ■) of a positive operator $ : H(di) — * Ti(d2) is a real-valued func- 
tion on the cone H+(d\), defined as the smallest concave function which coincides with 2\J a^ 2 ($(X)) 
on all rank 1 matrices X E H + (di). 

By concretizing the positive operator $ to be the partial trace, we can define these notions for 
bipartite matrices. The following definition is essentially from [5]. 

Definition 2.3. The concurrence C(M) of a positive semidefmite di ® d% bipartite matrix M is defined 

as C(tn; M), i.e. C(-) is the largest convex function on H+(did%) which coincides with 2\J a^ 2 (tri(X)) 
on all rank 1 matrices X E if+ (dis- 
similarly we introduce the following notion. 

Definition 2.4. The I-Gdelity F(M) of a positive semidcfinite d\ (8 di bipartite matrix M is defined 

as F(tn; M), i.e. F(-) is the largest convex function on H + (did2) which coincides with 2\J tr^ 2 {tr\ (X)) 
on all rank 1 matrices X E H+ididz). 

Since for any positive semidefmite bipartite rank 1 matrix ££* we have a^ 2 (tri(££*)) = av? 1 (^ r 2(CC*))j 
we arrive at the same definitions if the partial trace is taken with respect to the second subsystem [5] . 
Here £ E C dld2 denotes a vector. 

We now turn the the Lorentz cones and the associated Jordan structure. Let bq, ex, . . . , e„_i be the 
standard orthonormal basis vectors of R ra . For a vector x E M™, let Xq, x%, . . . , £„-i be the components 
of x with respect to this basis. Let further J„ be a diagonal matrix whose first diagonal element equals 
1 and whose all other diagonal elements equal — 1. 
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The the Lorentz cone L n c M n is defined as 



Ln = < X = (x 0) . . . ,X n -l) T G K™ | X > 



n-l 



\Y,4\- (2) 

\ fe=l 



The Lorentz cone can be viewed as the cone of squares for a certain Jordan algebra. From this Jordan 
algebra the spaces M" inherits the following structure [T] . 



Definition 2.5. For any vector x = (xq, . . . , x„_i) T € K n , let A± = xq± <J J2k=l x \ be the eigenvaJues 

of x. Define further A+ + A_ = 2xo to be the tracetr x and A+A_ = x^ — y~^l_^ x\ to be the determinant 
det x of x. 

With these definitions, a vector x is contained in L n if and only if its eigenvalues are nonnegative, 
and it is contained in the interior of L n if its eigenvalues are positive. Hence the Lorentz cone can be 
seen as an analogue of a 2 x 2 positive semidefinite matrix cone. The determinant of a vector x can be 
simply written as 

cist x — x JnX . 



We now define the linear map X : R 4 — > H(2) by 

T{x) = 



Xq +Xx X 2 + tX 3 
X 2 - ix 3 X - Xi 

As can easily be verified, the eigenvalues of the matrix I(x) equal the eigenvalues of x for any x £ R 4 , 
and T\L^\ = H + (2). Thus L4 can actually be identified with H+(2) by virtue of the isomorphism X. 

Remark 2.6. In the same manner, the cone L3 is isomorphic to the cone of positive semidefinite 2x2 real 
symmetric matrices and the cone L$ to the cone of positive semidefinite 2x2 quaternionic hermitian 
matrices. Therefore the formulae obtained in this paper equally apply for positive maps having as 
input space the space of real symmetric or quaternionic hermitian 2x2 matrices, given the definitions 
of concurrence and /-fidelity are adapted accordingly. 

We now extend the notions of concurrence and /-fidelity to M™- valued positive maps. Let E be a 
real vector space and K C E be a regular (i.e. closed, containing no lines) convex cone. The following 
notion is standard in convex analysis [8]. 

Definition 2.7. A point y £ K is said to be extremal if 1/1,1/2 £ K, y — yi + y 2 implies y\ = ctiy, 
y-2 = ctiy for some nonnegative scalars ct\, a 2 - 

The extremal points of the cone H + (d) of positive semidefinite matrices are precisely the matrices 
of the form where ( £ Cf The extremal points of the cone L n are precisely the points lying on 
the boundary dL n . 

Definition 2.8. Assume above notations and let m,n € N + be some integers. We call a linear map 
* : E -> R n K-to-L n positive if ^[K] C L n . We call a linear map * : R m -> R" Lorentz-positive if 
*[im] C L n . 

Definition 2.9. Let T be the set of extremal points of a regular convex cone K C E, where E is a 
real vector space. Let ^> : E — > K™ be a K-to-L n positive map, where n € N+ is some integer. The 
concurrence of ^ is the largest convex function C(^>; •) on the cone K which coincides with the function 
7(7) = 2-y/det(^(7)) on all extremal points 7 G T. 

Definition 2.10. Assume the notations of the previous definition. The I -fidelity of ^ is the smallest 
concave function F(^; ■) on the cone K which coincides with the function f(j) = 2-y/det($(7)) on all 
extremal points 7 G T. 

We thus define concurrence and /-fidelity as a convex and concave roof, respectively. In comparison 
to Definitions 12. 11 12~2"1 we do not restrict the input space of the operator ^ to be a space of self-adjoint 
operators, but in contrast we restrict its output space to be K™ equipped with a corresponding Jordan 
structure. This allows us to replace the second symmetric function a 2 by the determinant. As we will 
see further, this does in fact not at all restrict the generality. 

3 Relations between the different definitions 

In this section we investigate the relation between the different definitions of concurrence and /-fidelity. 
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First we explore the connection between the concurrence or /-fidelity of completely positive maps 
and that of bipartite matrices. 

Let $ : Tt(di) — > H(d 2 ) be a completely positive map, with Kraus representation 

${X) = Y,MXA* k . 

k=l 

With this representation we associate the third order tensor A Q ^ 7 , of dimension d± x d 2 X d% and with 
elements A Q ^ 7 = (Ay)^, a — 1, . . . , d±, = 1, . . . , d 2 , 7 = 1, ■ ■ • , c?3- Further we associate to this 
representation the d 2 d% x d\ matrix 

( A - \ 

A 2 



A 



V 4k J 



Then Q(X) equals the partial trace with respect to the first subsystem of the ds (g> d 2 bipartite matrix 
AX A*, $(X) = tn{AXA*). Therefore by Definitions O [U we have 

C($;X) = COijAXA*) = C(yLYA*). 

Now note that the concurrence of a bipartite matrix is independent of which subsystem is taken 
to define the partial trace, i.e. C(tri;AXA*) = C(tr 2 ; AX A*). If we define another positive map 
$' : H{d x ) -> H(d 3 ) by &(X) = tr 2 (AXA*), then it follows that 

C($;X) = C(&;X) 

for all X € H + (di). In a similar manner we obtain the equality 

F($;X) = F(&;X) V X e H + (d{) 

for the /-fidelity. 

The map has a Kraus representation given by 

&(X)=J^A' k XA' k *, (3) 
fe=i 

where the element (7,0;) of the d 3 x d\ matrix A'p is given by (Ay)p a = A a/ 3 7 . Hence the third 
order tensor associated to the Kraus representation ^ is given by ^ = A ai a 7 and has dimension 
d% x c? 3 x d 2 . The map $' can therefore be obtained from <E> by exchanging the last two indices in the 
corresponding third order tensor A. From the point of view of concurrence, the rank and the length 
of a completely positive map are hence interchangeable. We obtain the following result. 

Lemma 3.1. Let $ be a completely positive map with input space H.(di) 7 rank d 2 , and length d 3 . Then 
there exists a completely positive map $' with input space Ti(di), rank d%, and length d 2 such that 

C*($;X) = C(&;X), F($;X) =F(&;X) V X € H + {d x ). 

Moreover, for any matrix X G H + (di) of rank k < d\ there exists a d 2 ® d% bipartite matrix Mj of 
rank k such that 

C(fr,X) = C{M x ), F(§; X) — F(Mx)- 



The construction of $' and Mj is quite obvious from the above. This lemma may be a reason why 
the concurrence of completely positive maps of rank and length two |12j eventually boils down to an 
analogue of Wootters formula for the concurrence of 2 <g> 2 bipartite matrices [13] . 

We now explore the relation between Definitions 12.11 12.21 and 12.91 12.101 respectively. 
Let $ : H(di) — > H(d 2 ) be an arbitrary positive map. Let further {Mi, . . . , M d 2} be an orthonormal 
basis of the space Tt(d 2 ). We then have by (Q]) 

^($(X)) = M(tr$(X)) 2 -^($(X),M fc ) 2 ] . 
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Let us define a linear map $ L : W(di) -» R^ +1 by 

V2 



/ tr^(X) \ 



(4) 



Note that since $ is a positive map, we have $(X) >z and hence tr$(X) > ||$(X)|| 2 for all X > 0. 
Therefore the linear map takes the positive semidefmite matrix cone H+(di) to the Lorentz cone 
L d 2 +1 and is hence a H + (di)-to-L d 2 +1 positive map. With Definition ^. 51 we then have 

<j^(<S>(X)) = <P L (X) T J d 2 +1 <£ L (X) = det$ L (X). (5) 



The concurrences of the maps $, $l are defined by Definitions 12. 1[ 12.91 respectively. In view of the 
above equation they are the convex roofs of the same function. It follows that 

C{$;X) = C{$ L -X) VX&H+idx). 

In a similar manner we obtain 

F($;X) = F(<S> L ;X) VXeH+idi) 



for the /-fidelities, defined by Definitions 12 . 21 12 . 1 0[ respectively. We obtain the following result. 

Lemma 3.2. Let $ : 7i(di) — > 7i(d2) be a positive map. Then there exists a H + (di)-to-L d 2 +1 positive 
map <&l such that 

C($-X) = C{<S> L -X), F(<S>;X)=F(1> L ;X) V X e 

The map $l can be constructed as in Q). □ 
The lemma tells us that Definitions 12.91 12.101 are actually generalizations of Definitions 12.11 12.21 

4 Concurrence of Lorentz-positive maps 

There exists a simple case when we are able to compute the concurrence of a K-to-L n positive map 
explicitly, namely when the cone K in the input space is also a Lorentz cone. In this section we will 
derive a formula for the concurrence of Lorentz-positive maps. 

Let us first investigate the properties of certain real symmetric matrix pencils related to Lorentz- 
positive maps. 

Lemma 4.1. Let n > 3 be an integer and suppose that P — P T , J = J T are real symmetric n x n 
matrices, J is regular with signature (H — ■ • • — ) and there exists a number X € M such that P — A J >- 0. 
Let Ai, . . . , A„ be the generalized eigenvalues of the matrix pencil P — A J, with their real parts in 
decreasing order. Then the following assertions hold. 

(i) all eigenvalues Ai, . . . , A„ are real, 

(ii) Ai 7^ A2 and the matrix P — A J is positive definite if and only if X2 < A < Ai, 
(Hi) if A < Ai, then x T (P — \ J)x < implies x T Jx < for any x S K™. 

Proof. Since J is invertible, the generalized eigenvalues of the matrix pencil P — XJ are the eigenvalues 
of the matrix J _1 P and there are indeed exactly n of them. 

Since the signatures of J and — J are different, there must be at least one real number A such that 
P — XJ is singular, i.e. at least one generalized eigenvalue of the matrix pencil must be real. Denote the 
maximal real eigenvalue of J~ X P by A max and the minimal real eigenvalue by A m i n . Then for A > A max 
the matrix P — XJ has the same signature as — J, and for A < A m j n it has the same signature as J. 
It follows that the signatures of P — A J and P — XJ differ by exactly one sign for any A > A max and 
by exactly n — 1 signs for any A < A m ; n . Therefore the interval (A, A max ] contains at least 1 eigenvalue 
and the interval [A m i n , A) at least n — 1 eigenvalues. 

Thus all eigenvalues must be real, A max = Ai, A m i n = A„ and the interval (A, Ai) does not contain 
any eigenvalue. This proves (i), and (ii) follows by convexity of the cone of positive semidefmite 
matrices. 

Let us prove (iii). Assume the contrary, i.e. that there exists A < Ai and x G R n such that 
x T Jx > and x T (P - XJ)x < 0. Then we have x T {P - XiJ)x = x T {P - XJ)x + (A - X^x 7 ' Jx < 0. 
But P — Ai J >z by (ii), which leads to a contradiction and completes the proof. □ 
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We actually need the following version of this lemma with relaxed assumptions. 

Lemma 4.2. Let n > 3 be an integer and suppose that P = P T , J = J T are real symmetric n x n 
matrices, J is regular with signature (H — • • • — ) and there exists a number A £ R such that P — A J y 0. 
Let Ai, . . . , A„ be the generalized eigenvalues of the matrix pencil P — A J, with their real parts in 
decreasing order. Then the following assertions hold. 

(i) all eigenvalues Ai, . . . , A„ are real, 

(ii) the matrix P — A J is positive semidefinite if and only if A2 < A < Ai, 
(Hi) if A < Ai, then x T (P — XJ)x < implies x T Jx < for any x S R n . 

Proof. Assertions (i) and (ii) can be obtained by continuity arguments from the previous lemma when 
replacing P with P + eA for any e > 0, A y 0, and taking the limit e — > 0. Assertion (hi) is proven 
the same way as in the previous lemma. □ 

The following assertion is a well-known consequence of the 5-lemma [3"].|14|. 

Lemma 4.3. Let T : R m — > R™ be a Lorentz-positive map, represented by an n x m matrix. Then 
there exists A > such that such that T T J n T y XJ m . □ 
A proof is for instance in [4j Lemma 1]. 

This lemma relates the Lorentz-positive map T to the real symmetric matrix pencil T T J n T — A J m . 
By Lemma 14.21 this matrix pencil has m real eigenvalues. Let Ai, . . ., A m denote these eigenvalues in 
decreasing order. By assertion (ii) of Lemma \4. 21 A9 is the smallest A £ R such that T T J n T y XJ m . 
Hence by a convex separation argument there exists a nonzero x £ R m such that x T (T T J n T — A2 J m )x — 
and x T J m x < 0. As a consequence, x will be linearly independent from any vector in the interior of 
L m , since y T J m y > for any y e intL m . 

We shall now compute the concurrence of T, as given by Definition 12.91 Denote the positive 
semidefinite matrix T T J„T — A2 J m by Q. For any vector x € dL m we have x T J m x — and hence 



C*(T; x) = 2y/detT{x) = 2 Vx T T T J n Tx = 2^x T (T T J n T - A 2 J m )x = 2y / x T Qx. 

The concurrence C(T; •) : L m —>■ R is then the largest convex function on L m that coincides with the 
expression above on all vectors x € dL m , i.e. its convex roof (see [11]). Let us now show that the 
function 2\/x T Qx is actually this convex roof. 
For any x £ K m such that x T Qx > we have 



d 2 v / x T Qx _ Q{x T Qx) - Qxx T Q 
dx 2 ~ {x T Qxf/ 2 

If we evaluate this Hessian on the vector y, we get fa SgKg_Qg}ii^ Q g ) , Since Q y 0, we have 
(y T Qy)(x T Qx) > (y T Qx) 2 by the Cauchy-Schwarz inequality for the degenerate scalar product defined 
by Q and the Hessian is positive semidefinite. Thus the function 2y / 'x T Qx is convex in a convex 
neighbourhood U x of any vector x such that x T Qx > 0. However, if x T Qx = for some vector x, then 
the constant zero function supports 2y // x T Qx at x. It follows that the function 2y x T Qx is convex on 
the whole space M. m . 

Let us now prove that 2y/ ' x T Qx is a roof (for a definition see |llj). For any vector x £ intL m , 
define Lx as the 2-dimensional subspace L x C R m spanned by x, i. Since Qx = as a consequence 
of the relation x T Qx = and the positivity of Q, we have for any a, (3 £ R and x = ax + j3x that 
2\J x T Qx = 2\J u 2 x T Qx = 2\a\\/x T Qx. Hence the restriction of the function 2\J x T Qx on the set 
Sx = { x — ctx -I- (3x I a > 0} is linear. 

We shall now prove that S x H L m = L% n L m . Assume the contrary, i.e. that there exists a < 0, 
/3 £ R such that ax + (3x £ L m . Then (3 7^ 0, because ax £ — intL m . Moreover, —ax £ intL m and 
by convexity of L m we have h(— ax + (ax + (3x)) = ^x £ mtL m . But x 7^ ±intL m , which leads to a 
contradiction. 

Therefore any vector x £ int L rn can be represented as a convex combination of two points in 
dL m n S x . Hence the function 2\J x T Qx is indeed a roof. 

Note that x T T T J n Tx = det T(x) and x T J m x = det x. We have proven the following result. 

Theorem 4.4. Let T : R m — > R™ be a Lorentz-positive map. Then the pencil detT(a:) — A det a; of 
quadratic forms on R m has m real generalized eigenvalues. Let A2 be the second largest of them. Then 
the concurrence of T is given by the expression 

C(T;x) = 2v/detT(a;) - A 2 det a; 
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for all x £ L m . 

The function C(T;-) is linear along any affine subspace which is parallel to the eigenspace V 
corresponding to the eigenvalue A2 of the pencil detT(a;) — Xdetx. For any x 6 intL m there ex- 
ists a vector in V which is linearly independent of x. This allows to obtain an optimal decomposi- 
tion of x as convex combination x = fiy + (1 — fi)z of two points y,z £ dL m such that C(T;x) = 
/iC(T;y) + (l-/i)C(T;z). ' □ 



5 /-fidelity of Lorentz-positive maps 

In this subsection we prove a result analogous to Theorem 14.41 for the /-fidelity of Lorentz-positive 
maps. Assume the notations of the previous section. 

Theorem 5.1. Let T : W n — > R™ be a Lorentz-positive map. Then the pencil detT(a;) — A det 2 of 
quadratic forms on R m has m real generalized eigenvalues. Let A m ; n be the smallest of them. Then the 
I -fidelity of T is given by the expression 

F(T; x) = 2^detT(a;) - A min det 2; 

for all x £ L m . 

The function F(T; •) is linear along any affine subspace which is parallel to the eigenspace V corre- 
sponding to the eigenvalue A m i n of the pencil det T(x) — A det x. For any x £ int L m there exists a vector 
in V which is linearly independent of x. This allows to obtain an optimal decomposition of x as convex 
combination x = fj/y+(X — /j,)z of two points y, z £ dL m such that F(Y; x) = [j,F(T; y) + (1 — fi)F(T; z). 

Proof. Denote the generalized eigenvalues of the matrix pencil T T J„T — \ J m by Ai, . . . , A m , in de- 
creasing order. By Lemma 14.31 Lemma 14.21 is applicable and all generalized eigenvalues of the pencil 
are real. 

Define the matrices Q m = T T J„T — X m J m and Q\ — T T J„T — Ai J m . The /-fidelity of T at a 
vector x £ dL m is given by 

F(T; x) = 2v/detT(x) = 2y/x T T T J n Tx = 2^Jx T (T T J n T - X m J m )x = 2y / x T Q m x, 

because x T J m x = for any x £ dL m . The /-fidelity F(T; •) : L m — > K is then the smallest concave 
function on L m that coincides with the expression above on all vectors x £ dL m , i.e. its concave roof 
(see [H]). We have to show that the function 2\J x T Q m x is actually this concave roof. 

We have Q m = Q± + (Ai — A m ) J m . Hence for any vector x £ L m we have x T Q m x = x T Q\x + (Ai — 
X m )x T J m x > 0, because Q\ > by assertion (ii) of Lemma [4~2l Ai > A m and x T J m x > 0. Hence 
2y // x T Q m x is well-defined on L m . 

If Ai = A m , then all generalized eigenvalues of the pencil are equal and the polynomial p(X) = 
det(T T J„T — XJ m ) has an m-fold root at A = X m . Hence the matrix Q m must be zero, in which case 
2^J x T Q m x = is a concave roof. 

Let us henceforth assume that Ai > A m . Let x be an eigenvector to the eigenvalue X m of the pencil 
T T J„T — XJ m . Then Q m x = 0, x T Q m x = and x T J m x < by assertion (iii) of Lemma l4.2l It follows 
that x^x are linearly independent for any x £ int/v m . 

Let us show that 2\J x T Q m x is concave on L m . For any e > the signature of the matrix Q m + £J m 
equals that of J m , because Q m + eJ rn is regular for any e > 0. Consequently, Q m has at most one 
positive eigenvalue and we can represent Q m as S T J m S, where S is some singular real m x m matrix. 
Let now x £ intL m . Then we have x T J m x > and x T Q m x = x T Q\x + (Ax — X m )x T J m x > 0. As in 
the previous section, it follows that 

d C^jyiX S Jm S ' x S J jyiS x S JfuSxx S Jm S 
dx~ 2 = (x T Q m x) 3 / 2 ' 

Let us evaluate this Hessian on the vector y. Denote Sx by x and J m Sy by y. Then we get 

T d 2 y / x T Q m x _ y T J m y ■ x T J m x - y T x~x T y _ y T (xx T - x T J m x ■ J m )y 
V dx 2 V ~ {x T Q m xf/ 2 ~ {x T Q m xf/ 2 

Note that x T J m x = x T Q m x > and hence x £ ±intL m . It follows that J m 5 £ ±inti m . But J m x is 
an eigenvector of the quadratic form xx T — x T J m x ■ J m with eigenvalue zero. 

Assume that this form has an eigenvector z with negative eigenvalue. Then z is orthogonal to J m x 
and hence z T J m z < 0. This is because a vector in intL m cannot be orthogonal to any other nonzero 
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vector in L m . But then z T (xx T — x T J m x ■ J m )z = (x T z) 2 + (x T J m x)(— z T J m z) > 0, which leads to a 
contradiction. 

Thus xx T — x T J m x ■ J m >z and - ^q x ® — ^ 0. It follows that the function 2\J x T Q m x is concave 
on intL m and hence by continuity on the whole cone L m . 

By the same arguments as for the proof of Theorem 14.41 any vector i 6 int L m can be represented 
as a convex combination of two points in dL m fl Ij = dL m n Sx, where L x , Sx are defined as in the 
previous section and the restriction of the function 1\J x T Q m x on the set Sx is linear. Thus 2a/ x T Q m x 
is indeed a roof. □ 



6 Concurrence and /-fidelity of positive maps 

In this section we apply the results of the previous sections to compute the concurrence and the 
/-fidelity of positive maps with input space Ti.(2) and bipartite density matrices of rank 2. 

As we have seen in Section 2, the cone H + (2) of positive semidefinite complex hermitian 2x2 
matrices is isomorphic to L4, with the matrix determinant being equal to the determinant induced on 
M 4 by the Jordan structure of L4. By application of Lemma EOl and relation Theorems 14.41 and 15.11 
then yield the following results. 

Theorem 6.1. Let $ : Tt(2) — > 7i(d 2 ) be a positive operator. Then the pencil ^(^{X)) — AdetX of 
quadratic forms on TL{2) has 4 real generalized eigenvalues. Denote these eigenvalues by Ai, A2, A3, A4, 
in decreasing order. Then the concurrence and the I -fidelity of 4> are given by the expressions 

C($; X) = 2a/^ 2 ($(X)) - A 2 detAT, F(<$>; X) = 2y/a% 2 ($(X)) - A 4 detX 
for allX e H+(2). 

The function C(<I>; •) (respectively, -F( C E > ; ■)) is linear along any affine subspace which is parallel to 
the eigenspace V2 (respectively, V4) corresponding to the eigenvalue A2 (respectively, A4J of the pencil 
o~2 2 (<&(X)) — AdctA. For any X G intiJ + (2) there exists a matrix in V% (respectively, V4) which is 
linearly independent of X . This allows to obtain an optimal decomposition of X as convex combination 
X = nY+(l-fi)Z of two rank 1 matrices Y, Z <E H + {2) such thatC(<£;X) = nC($; y) + (l-,u)C($; Z) 
(respectively, F(<f>; X) = Y) + (1 - Li)F(<f>; Z) ). □ 

Note that the formulae provided by Theorems 14.41 and l 5 . 1 1 are coordinate independent. Moreover, the 
formulae do not change if the quadratic form det x is scaled by multiplication with a positive number. 
We are therefore not bound to the standard Lorentz cone L n as defined in ((2]). All we need are two 
quadratic forms on the input space, one having signature (H — • • • — ) and defining a convex cone by its 
zero set, and the other describing the value of the concurrence or /-fidelity on the extremal elements of 
this cone. In fact, the proofs of Theorems 14.41 and 1 5 . 1 1 used the assumption that T is Lorentz-positive 
only to ensure the assumptions of Lemma 14.21 We can hence relax the Lorentz-positivity of T and 
even abstract ourselves from the notion of a positive map. We have the following general result, whose 
proof goes along the lines of proof of Theorems 14.41 and 1 5 . 1[ with obvious modifications. 

Theorem 6.2. Let E be a real vector space of dimension n and let J be a regular quadratic form on 
E with signature (H — •••—). The set {x G E \ J(x) > 0} forms two convex cones, which are linearly 
isomorphic to the Lorentz cone L n . Let K C E be one of these two cones. Let P be another quadratic 
form on E satisfying the condition 

3 A e E : P h XJ. (6) 

Then the pencil P — XJ of quadratic forms has n real generalized eigenvalues. Let X±, . . . , A„ denote 
these eigenvalues in decreasing order. Let the function p : 5K — > K be defined by p(x) = 2yJP{x). 

Then the largest convex function p CO nv and the smallest concave function p conc on K which coincide 
with p on dK are given by 

Pconv(x) = 2 V /P(.t) - A 2 J(ir), Pconc{x) = 2y/P(x) - X n J(x). 

These two functions are hence the convex and the concave roof of p, respectively. 

Let V2, V n C E be the eigenspaces corresponding to the eigenvalues A 2 , X n , respectively, of the pencil 
P — A J. Then p con v is linear on any affine subspace parallel to V 2 and p conc is linear on any affine 
subspace parallel to V n . For any x £ int K there exist elements in V2 and V n which are linearly 
independent of x. This allows to obtain optimal decompositions of x as convex combinations x = 
/iyi + (l-^)y 2 = f]Zi + (l-r])z2 of points yi,y 2 ,zi,z 2 € dK such that p conv (x) = /ip(yi) + (1 - fi)p(y 2 ), 
Pconc(x) — r}p(zi) + (1 - rj)p(z 2 ). □ 
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Let di > 2 and let P C C dl be a linear complex subspace of dimension 2. Denote by Up the 
subspace of all matrices in TL{d\) whose range is contained in P. Then Up is isomorphic to 7i(2) and 
the intersection Kp = Up D H + (d\) is isomorphic to H + {2) and hence to L4. Moreover, if an element 
of Kp is represented as convex combination of extremal elements of H+(di), then all these extremal 
elements also have to lie in Kp. Therefore the concurrence of an element X E Kp with respect to 
a positive operator <E> : H(di) — > H{d 2 ) equals the concurrence of X with respect to the restriction 
of $ to Up. The same holds for the /-fidelity. Thus Theorem 16.11 is applicable also in this case, but 
we have first to find a function on Up which can serve as determinant in the sense of Definition 12.51 

It suffices that this function is a quadratic form with signature (H ), is positive on int -fTp and 

zero on dKp. The second symmetric function 

CT2 1 on TC(di) has signature (H — ■ • • — ), is zero on all 
extremal elements of H + (di) and positive on all other elements of this cone. Hence its restriction to 
Up fulfills the necessary requirements. We can deduce the following result. 

Theorem 6.3. Let $ : Tt(di) — > TL{d2) be a positive operator and let X € H + (di) be a matrix of rank 
not exceeding 2. Let further P C C dl be a linear complex subspace of dimension 2 such that the range 
of X is contained in P. Denote by Up the subspace of all matrices in TL{d{) whose range is contained 
in P. 

Then the generalized eigenvalues of the pencil a d2 o <&\u P — AcTj \u p are a ^ real. Denote them by 
Ai,A2,A3,A 4 in decreasing order. Then the concurrence and the I-fidelity of X with respect to $ are 
given by 

C(*; X) = 2^/a d 2 2 (<i>(X))-X 2 a d 2 1 (X), X) = 2^ (4>(X)) - X^ (X). 

The function C($; ■) (respectively, F(<f>;-)) is linear along any affine subspace which is parallel to 
the eigenspace V2 ( respectively, V4 ) corresponding to the eigenvalue A2 ( respectively, A4 ) of the pencil 
a 2 2 o <5>\u P — Xa 2 x \u P - If X has rank two, then there exists a matrix in V 2 (respectively, V4) which is 
linearly independent of X . This allows to obtain an optimal decomposition of X as convex combination 
X = nY+{\-n)Z of two rank 1 matrices Y,Z € H + (di) such that C($; X) = fiC(§; Y)+(l-n)C(®; Z) 
(respectively, X) = ^F($; Y) + (1 - ii)F($; Z) ). □ 

For a d\ ® d 2 bipartite matrix X the concurrence is defined as concurrence of X with respect to 
the partial trace. Setting $ to the partial trace in the previous theorem, we obtain a formula for the 
concurrence of bipartite matrices. 

Corollary 6.4. Let X be a d\ ® d 2 bipartite matrix of rank not exceeding 2. Let further P C <C dld2 be 
a linear complex subspace of dimension 2 such that the range of X is contained in P. Denote by Up 
the subspace of all matrices in Ti(did2) whose range is contained in P. Define the two quadratic forms 
Q-L-.Ah^ 2((fr A) 2 ~ tr((triA) 2 )) and Q 2 : A^ (tr A) 2 - tr A 2 on H(did 2 ). 

Then the generalized eigenvalues of the pencil Qi\u P ~~ ^Q2\u P are oil real. Denote them by 
Ai, A2, A3,A 4 in decreasing order. Then the concurrence and the I-fidelity of X are given by 

C(X) = ^Qi(X)-\ 2 Q 2 {X), F(X) = VQ1PO-A4Q2PO. 

The function C($; ■) (respectively, F(Jf>;-)) is linear along any affine subspace which is parallel to 
the eigenspace V 2 (respectively, V<t) corresponding to the eigenvalue X 2 (respectively, A4J of the pencil 
Qi\u P ~ XQ 2 \u P . If X has rank two, then there exists a matrix in V 2 (respectively, V4) which is 
linearly independent of X . This allows to obtain an optimal decomposition of X as convex combination 
X = /jY+(1- [x)Z of two rank 1 matrices Y, Z G H + (dxd 2 ) such that C(X) = fJ,C(Y) + (1 - fJ,)C(Z) 
(respectively, F(X) = (iF(Y) + (1 - fi)F(Z) ). □ 

The quadratic form Q\ in the corollary can be replaced by the form A 1— > 2((tr A) 2 — tr((tr 2 A) 2 )) 
or by ® Sd 2 , where Sd is the universal inverter on Tt(d) as defined in [5]. All these forms are equal 
on the pure states and hence on the boundary of Kp. 

Remark 6.5. While the form A 1— > 2((tr A) 2 — tr((tr 2 A) 2 )) is generated by the determinant of the 
Lorentz-positive map tr 2 and hence fulfills condition ([6]), the argument for ® 5<j 2 , or for any other 
quadratic form that coincides with Qi on the boundary of Kp, is as follows. Let J be a regular quadratic 
form on W 1 with signature (H — • • ■ — ). Let further Q, Q' be quadratic forms fulfilling x T Qx = x T Q'x 
for all x such that x T Jx = 0. Then there exists a real number A such that Q' = Q + XJ. Therefore 
Q — XJ and Q' — XJ are essentially the same pencils. 

7 Density matrices of 2-edge graphs 

In this section we apply the results of the previous section to density matrices of graphs with 2 edges. 
Let G be a graph and M(G) its adjacency matrix, i.e. a symmetric matrix of size n x n, where n is the 
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number of vertices of G, its {k, I) element given by 1 if the vertices k, I are connected by an edge and 
by otherwise. Let further A(G) be the degree matrix of the graph, i.e. the unique diagonal matrix 
such that the combinatorial Laplacian L(G) = A(G) — M(G) contains the vector (1, 1, ... , 1) T in its 
kernel. The density matrix p G of G is defined as tr L ^ = 2n , where n e is the number of edges. For 
details see [2] . Note that n e is an upper bound for the rank of the density matrix. 

If the vertices of the graph are arranged in an array of size d\ x c?2, then its density matrix inherits 
a natural d\ ® cfe bipartite structure. In [5] the formula of Wootters [13] was used to compute the 
concurrences of all graphs on 4 vertices, arranged in an 2 x 2 array. It turned out that either the density 
matrix was separable or the concurrence was equal to The corresponding optimal decomposition 
of the nonseparable density matrices derives from the decomposition of the combinatorial Laplacian in 
the combinatorial Laplacians of all subgraphs with one edge. 

In this section we compute the concurrences and /-fidelities of all density matrices of graphs having 
rank 2. These are the density matrices of all graphs with 2 edges and of all graphs with 3 edges arranged 
in a closed loop. For the computation we used the universal inverter as defined in [S]. Suppose the 
vertices of the graph G in question are arranged in an array of size d\ x di , giving rise to a d\ <g> di 
bipartitioned density matrix pc- Let the form Q\ on "Hidydz) from Corollarv l6.4l be defined by S& x (&Sd 2 , 
i.e. acting as Qi : pa *-* (pg, tr pGldKh ~ Id 1 ®tr 1 p G ~tr2PG®Id 2 +Pg), and let the form Q 2 be defined 
by Q2 ■ Pg !— ► {t r Pg) 2 — tr Pg- Let further U be the 4-dimensional subspace of 7i{d\di) consisting 
of all matrices that have the same range as pc- Then the concurrence and the /-fidelity of pc are 
given by C(p G ) = \/Qi{pg) ~ A 2 Q 2 (pg), F(p G ) = ^Qi{pc) - MQ2{pg), where Ai, A 2 , A 3 , A 4 are the 
generalized eigenvalues of the pencil Qi\u — \Q2\u in decreasing order. 

Below we depict the graphs, each representing one isomorphism class, and list the values of the 
corresponding concurrences and /-fidelities, the generalized eigenvalues of the pencil <5i |c/ — AQ2 1 c/ and 
the values of the quadratic forms Qi, Q2 on pc- 
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The graphs of subtypes distinguished by letters lead to density matrices which have the same range 
and hence share the space U and the spectrum Ai, . . . , A4 of the corresponding pencil. 

8 Conclusions 

The concurrence of positive operators or bipartite positive semidefinite matrices is defined as the 
convex roof of the square root of a certain quadratic function defined on the positive semidefinite rank 
1 matrices in the input space TL{d), the space of complex hermitian matrices of size d x d. The convex 
roof is the largest convex extension of the function in question to the convex hull of the rank 1 matrices, 
i.e. to the cone of positive semidefinite matrices in Tl(d). Our ability to compute concurrences hence 
depends on our ability to compute such convex roofs. 

The concept of convex roof can be generalized to other regular convex cones than the positive 
semidefinite cone. It can be defined as the largest convex extension of a function defined on the 
extremal elements of the cone to the cone itself. In this paper we derived an explicit expression for the 
convex roof of a certain class of functions and a certain class of cones, namely cones that are generated 
by quadratic forms with signature (H — •••—), i.e. linear images of the Lorentz cone. For the same 
class we are also able to compute the concave roof. This result is the main technical contribution of 
the present paper and is formalized in Theorem 16.21 

It allows us to obtain explicit formulae for the concurrence of Lorentz-positive maps as defined 
in Definition 12.91 as well as for a function called /-fidelity, which is defined via the concave roof in 
Definition 12.101 These formulae are provided in Theorems 14.41 and 15.11 From Lemma 13.21 it follows 
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that if the input space of a positive operator is 7i(2), then computing its concurrence and /-fidelity 
can be reduced to computing the concurrence or /-fidelity of a Lorentz-positive map. This allows us 
to obtain explicit formulae for these quantities, provided in Theorem 16. II More generally, it allows us 
to compute the concurrence and /-fidelity of any matrix of rank not exceeding two with respect to any 
positive operator. These formulae are provided by Theorem 16.31 Further we investigate the relation 
between the concurrence of positive operators and bipartite matrices. In Lemma 13.11 we show that 
the notions of concurrence and /-fidelity for a bipartite matrix are essentially equivalent to those for 
a completely positive operator. Hence we are also able to compute the concurrence and /-fidelity of 
rank two bipartite matrices. The corresponding formulae are provided in Corollary 16.41 In all cases 
the optimal decomposition yielding the value of the concurrence or /-fidelity contains two pure states. 
The optimal decomposition can be obtained via the eigenvector of the corresponding matrix pencil to 
the second largest generalized eigenvalue for the concurrence and the smallest one for the /-fidelity. 
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